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ABSTRACT 

The origin of the Galactic center diffuse X-ray emission (GCDX) is still under intense inves- 
tigation. In particular, the interpretation of the hot {kT « 7keV) component of the GCDX, 
characterised by the strong Fe 6.7 keV line emission, has been contentious. If the hot component 
originates from a truly diffuse interstellar plasma, not a collection of unresolved point sources, 
such plasma cannot be gravitationally bound, and its regeneration would require a huge amount 
of energy. Here we show that the spatial distribution of the GCDX does not correlate with the 
number density distribution of an old stellar population traced by near-infrared light, strongly 
suggesting a significant contribution of the diffuse interstellar plasma. Contributions of the old 
stellar population to the GCDX are implied to be ~ 50 % and ~ 20 % in the Nuclear stellar disk 
and Nuclear star cluster, respectively. For the Nuclear stellar disk, a scale height of 0?32 ± 0?02 
is obtained for the first time from the stellar number density profiles. We also show the results of 
the extended near-infrared polarimetric observations in the central 3° x 2° region of our Galaxy, 
and confirm that the GCDX region is permeated by a large scale, toroidal magnetic field as previ- 
ously claimed. Together with observed magnetic field strengths close to energy equipartition, the 
hot plasma could be magnetically confined, reducing the amount of energy required to sustain it. 

Subject headings: Galaxy: center — X-rays: ISM — polarization — ISM: magnetic fields 



1. INTRODUCTION 

In the late 80s, a Japanese X-ray satellite 
GINGA revealed the presence of a diffuse and 
rather uniform 6.7 keV emission from highly ion- 
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ized. Helium-like ions o f iron at the direction of the 
Galactic center (GC; iKovama et al. 19891 ). The 
line emission and associated continuum compo- 
nent, called the Galactic center diffuse X-ray emis- 
sion (GCDX), resembles the G alactic ridge diffuse 



X-ray emission (GRXE; e.g., ICooke et all Il969l: 



Kovama et al.l 119861 ) extending more than 100° 
along the Galactic plane. For the GRXE, more 
than 80 % of the diffuse emission has been claimed 
to be resolved into point sources teevnivtsev et al" 



20091 ). suggesting faint X-ray point sources in ori- 



gin. 

On the other hand, the origin of the GCDX, in 
particular its very hot component with a temper- 
ature of kT ^ 7keV, is more puzzling. So far only 
10 - 40 % of the GCDX has been claimed to be re- 
solved into fa int point sources even with the Chan- 
dra s atellite ( Muno et al. 2004; Revnivtsev et al.l 
20071 ). The plasma temperature, represented by 



1 



flux ratios of iron emission lines, is systemati- 
cally higher for the G CDX than for the GRXE 
(|Yamauchi et al.ll2009l ). indicating their different 
origins. 

Two main ideas have been suggested to account 
for it : a truly dif fuse plasma tha.t bath es the emit- 
ting region (e.g., Kovama et al.l [l989l ) : and a su- 
perposition of a large number of unresolved poin t 
sources as the GRXE is fe.g.. IWang et al.l 120021) . 
In the later case, candidates are old stellar bi- 
nary systems such as cataclysmic varia bles (CVs 



and c oronally active binaries (ABs; Sazonov et al 



20061 ). So if the hot component originates in the 
discrete sources, its spatial distribution should be 
very similar to that derived by old stars observable 
in infrared wavelengths. For this purpose, a stellar 
mass distribution model constructed from infrared 



surface bright ness maps (Launhardt et al.M2002r) 
has been used iMuno et aL 20091 : Uchivama et al 



20111: iHeard fc Warwickl l2013l) However, such 



maps could be subject to the influence of bright 
stars. The angular resolution in these maps was 
only 0?7, so that the stellar density profile in the 
direction orthogonal to the Galactic plane had 
to be inferred from proxies (dust emission, radio 
emission from molecular clouds). In addition, un- 
certainties of the mass model seem to be as high 
as a factor of two ( Launhardt et al. I I2OO2I) . 

We have constructed a stellar number density 
map of the GC region from new near-infrared 
(NIR) observations with more than 1,000 times 
higher spatial resolution (Yasui et al. in prepara- 
tion), which enables us to directly compare the 
stellar distribution with GCDX. In this letter, 
we summarize the NIR imaging observations and 
their results, and provide additional evidence for 
the hypothesis that the GCDX arises from a truly 
diffuse hot plasma. We also show results of our 
recent polarimetric observations. The results pro- 
vide strong evidence for a large-scale toroidal mag- 
netic field configuration which could confine the 
hot plasma magnetically. 



the 1.4 m telescope IRSF. SIRIUS provides J (1.25 
/im), H (1.63 /xm), and Ks (2.14 /im) images si- 
multaneously. The averages of the lOcr limiting 
magnitudes are H = 16.6 and Ks — 15.6. We 
do not use the J-band data due to severe inter- 
stellar extinction. Fu rther details are given in 
Nishivama et al.l (|2006[ ). 

The stellar number density map is constructed 
as follows: At first, an H and Ks color magni- 
tude diagram (CMD) is constructed for each sub- 
field of 20' X 20', and foreground sources with their 
blue H — Ks color are removed. The typical color 
of stars in the GC is H — Ks > 1.0, and color 
cuts to remove the foreground sources are 0.3 - 
1.1. We carry out an extinction correction for 
each star using the observed H — Ks color, the 
mean intrinsic color of {H — Ks)o « 0.20 (con- 
sidering t he limiting magnitude s and the Galactic 
model bv lWainscoat et a~ 1992 ). and an interstel- 
lar extinction la w, A(Ks) = 1 .44 x {{H - Ks) - 
{H — Ks)o} (Nishivama et al.l i2006). Here we ob- 
tain an extinction-corrected Ks-hand magnitude, 
Ks,o, and the amount of interstellar extinction for 
each star. 
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Fig. 1. — The stellar number density map of the 
central 6° x 2° region of our Galaxy. Stars with 
Ks,o < 8.0 and < 10.5 are used for the central 20', 
and for the outside of 20', respectively. Several 
low-density regions are seen in the NB, and are 
not used in the following analysis. 



2. Observations and Data Analysis 

The central region of our Galaxy, | I |< 3?0 
and I b |< 1?0 (corresponding to 840 pc x 
280 pc at 8kpc from the Sun), was observed 
from 2002 to 2004 using the NIR camera SIRIUS 
( Nagashima et al. 1999t Nagavama et al on 



We then construct a stellar number density 
map using stars with i^s.o < 10.5 (Fig. [1]). 
For the central 20' x 20' field, source confusion 
is so severe that a different magnitude limit of 
Ks,o < 8.0, and a conversion factor derived by 
the ratio of the number of stars with Ks^o < 8.0 
and Ksfi < 10.5 are used. By comparing Ks- 



2 



band luminosity functions constructed with the 
extinction-corrected sta rs and a Galactic model 
(jWainscoat et alJ Il992l ) , these magnitude limits 



are determined so that completeness at the limits 
is almost 100 %. We make completeness correc- 
tions with recovery rates determined by adding ar- 
tificial Ks = 12.5 sources in the images, and con- 
firm that th e recovery rates are very high, ~ 96 % 
in average ( Hatano et al.l l2013[ ) . For more de- 
tail, see Yasui et al. (in preparation). Also, at 
the very center (< 1'), we have used images ob- 
tained with an 8-m telesco pe VLT and ISAAC 
(|Nishivama fc Schodelllioii . 

To determine a large-scale interstellar MF con- 
figuration, we have carried out NIR polarimetric 
observations usin g IRSF and a NIR polarimetric 
imager SIRPOL |Kandori et al.ll2006l ). from 2006 
to 2010. We have extende d the survey region from 



I I \< 1?0 (and I b |< 1?0: iNishivama et al.l[201ol ) 
to I I 1^ 1°5, which covers almost the whole region 
dominated by GCDX. Comparing the polarization 
between stars distributed further and closer side in 
the GC, we obtain polarization originating from 
magnetically al igned dust grains in the G C (for 



more detail, see lNishivama et al.|[2009ll2010l) . The 

polarized angle traces the GC's MF direction pro- 
jected onto the sky. 

3. Results and Discussion 

3.1. Stellar Number Density Profiles 

A disk-like structure is seen in the stellar num- 
ber density map (Fig. [T]). This is known as the 
Nuclear bulge (NB), which consists of the Nuclear 
stellar disk (NS D) and the Nucl ear stellar clus- 



ter (NSC; iLaunhardt et al.l 120021 ). Here our ob- 



servations clearly reveal morphology of the NB 
on a large scale, with much higher spatial reso- 
lution than previous studies. The NB has a sym- 
metric, disk-like structure with a scale height of 
0?32 ± 0?02 (Fig. HI), although several low-density 
regions are seen. Those are very dense molecu- 
lar clouds in front of /inside the NB. A true stellar 
number density is difficult to be derived in these 
regions, and thus they are not used in the following 
analysis (see Yasui et al. in preparation). 

The longitudinal and latitudinal profiles of 
the 6.7 keV line emission measured by Suzaku 



(jKovama et al.ll2007l : IUchivama et al.ll201l[ ) clearly 
show an excess at the NB region over the stellar 
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Fig. 2. — Top: Longitudinal profile of the stellar 
density distribution after a completeness correc- 
tion. Only data points at < is plotted, and 
those with pink marks are not used for the fit- 
ting. The profile can be approximated by a sum 
(black line) of two different power laws (blue and 
cyan lines) and a constant component (red line). 
The boundary of the two power laws, /* = 0?63, 
corresponds to the out er edge of the inner NB 



(jLaunhardt et al.l 120021 ). Bottom: Latitudinal 
stellar density profiles. Data points from posi- 
tive and negative 6* are plotted as green (cyan) 
and red (pink) marks, respectively. The cyan and 
pink marks are not used for the fitting because the 
number densities are underestimated in the corre- 
sponding regions due to strong line-of-sight extinc- 
tion. Data points at | 6, |< 0?1 are also not used 
for fitting because of the dominance of the NSC. 
The profile is fit with the NSD (blue line) and the 
Galactic disk (GD; red line) components using ex- 
ponential functions, where the scale height of the 
GD is assumed to be 1?3. The NSD component 
has a scale height of 0?32 ± 0?02. The number 
densities are calculated in rectangular bins with a 
size of l'{l) X 2' (6) for the longitudinal profile at 
6* = 0°, and 2'{l) x l'(6) for the latitudinal profile 
at h = 0°. 
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mimber density profiles (Fig. [3]). The two profiles 
are overplotted and scaled to have the same values 
at 1?5 <| /* |< 2?8, i.e., in a region outside of the 
NB [U and 6* denote the angular distance from 
Sgr A* along the Galactic longitude and latitude, 
respectively, and 6*) = (/ + 0?056, b + 0?046)]. 
Fitting the longitudinal profile in the range 
-0?7 < < -0?1 with a power-law of oc 6*"", 
where is angular offset from Sgr A*, gives 
Q^star = 0.30 ± 0.03 for the stellar number den- 
sity. This is different both from 0.44 ± 0.02 for 
the 6.7 keV profile in the same range, and from 
0.60^0 03 for the integrated emission o f Fe6.7 and 



6.9 keV lines (|Heard fc Warwickll2013l ). 

The majority of faint X-ray (L2-iokcV < 
1030 ergs-i) sources which have not been re- 
solved but contribute to the GCDX are most 



likely to be old binary systems (|Sazonov et al 
2006). Using the synthe tic CMD computation 
(jAparicio fc Gallarti |200J) , and a constant star 
formation history during 13Gyr for the NSD 



KV 

(jFiger et al.ll2004l ). we have confirmed that about 



75% of the stars with i^s.o < 10.5 are older than 
1 Gyr. So the NIR map and profiles shown here 
trace the distribution of the old stars, and they are 
clearly different from those of the 6.7 keV emission. 
A contribution from faint discrete sources to 



the GCDX has been claimed (IWang et al.M2002 : 
Muno et al.l 120041 : iRevnivtsev et al.ll2007h . To in- 
vestigate the contribution of the point sources, es- 
pecially of the old stellar population detectable in 
our observations, we construct longitudinal and 
latitudinal profiles for the ratio of the 6.7 keV 
emission to the stellar number density (Fig. 
When the profiles are scaled to be unity at 1?5 <| 
h \< 2?8, the ratios are 1.5 and - 3 in the NSD 
and NSC, respectively. 

The NB and NSC have a different forma- 
tion history from the Galactic bulge (GB), and 
have formed stars over their entire lifetime, in- 
dicating more bright stars in the NB and NSC. 
Here we use the synthe tic CMD computation 
( Aparicio &: Galla"r3l2004 ) to estimate the fraction 
of the total number of stars formed in the compu- 
tation, TVaii, to the number of stars with K < 10.5, 
Nk<io.5- This ratio, R = Ni,ii/Nk<io.5, repre- 
sents the ratio of the theoretically expected total 
number of stars to the number of stars detected 
in our observations. The star formation histo- 
ries used here are: a burst star formation from 
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Fig. 3. — Longitudinal (top) and latitudinal (bot- 
tom) profiles of the stellar number density after 
a completeness correction (red crosses). Over- 
plotted are the 6.7keV-emission profiles (b lue x; 
Kovama et"aLll2007t lUchivama et allboill ). The 
region outside the NB, 1°5 <\ h \< 2?8, is used 
to scale the the 6.7 keV emission profile to have 
the same value as the stellar number density. The 
number density is c alculated in the same r ectan- 
gles as those used in lUchivama et ah with 
a size of 0?1(Z) x 0?2(6) for the longitudinal profile, 
and 0?2(0 x 0?1(6) at the position of I = -0?17 
for the latitudinal profile. The same scaling factor 
is used for the latitudinal profile. 
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Fig. 4. — Longitudinal (top) and latitudinal (bot- 
tom) profiles of the ratios of the 6.7 keV emission 
to the stellar number density, scaled to be unity 
at the position for normalization, 1?5 <| |< 2? 8. 
These profiles represent a contribution of point 
sources, traced by our NIR observations, to the 
GCDX in the assumption that the contribution of 
truly diffuse hot plasma is negligible at the po- 
sition for normalization (i.e., the Galactic ridge 
region) . 



10 to 13 Gyr ago fo r the GB (i.e., outside the NB; 
Zoccali et al 12003 1): a constant star formation rate 
for 13 Gyr for the NSD (Figcr et al. 2004); and the 
history derived bv [Pfuhl et al.l (1201 ll) for the NSC. 
When we scale the ratio R to be unity for the GB, 
we obtain i?GB : ^^nsd : ^nsc ~ 1 '■ 0.8 : 0.6. 
This result means that the densities of the old 
stellar population in the NSD and NSC are over- 
predicted by the bright stars. Taking into account 
this ratio, we have found that the contributions 
of the old stellar population to the GCDX are 
(1/1.5) x 0.8 ~ 0.5 and (1/3) x 0.6 - 0.2 for 
the NSD and NSC, respectively, in the assump- 
tion that the X-ray luminosity function is uni- 
versal, and that the contribution of point sources 
to the GRXE is 100%. The contribution at the 
NSC is in goo d agreement with ~ 1/6 derived by 
Kovama et all ((2009). 

A larger X-ray emissivity per unit stellar 
mass for the GCDX than the GRXE has been 
claimed to explain the di fferent distributions of 
stars and 6.7 keV emissio n ( Revnivtsev et al.f2007 ; 
Heard fc Warwick! 12013'). To change the emissiv- 
ity, at least one of initial mass function (IMF), 
binary fraction (BF), or star formation history 
(SFH) is required to be different in the NB from 
the GB. In the preceding paragraphs, we have 
shown that different SFHs cannot explain the dif- 
ferent spatial distributions of the stellar number 
density and 6.7 keV emission. Co nsidering an uni- 
versal IMF ( Bastian et al. 2010fl and a po ssible 
top-heavy IMF in the GC (|Figer et al.lll999t ). the 
number of old, low- mass stars (i.e., CVs and ABs) 
per unit stellar mass never increases, it only de- 
creases. Also, a higher stellar density tends to 
destroy binaries rather than to form them via a 
capture process, whi ch seems to pl ay a small role 
in binary formation (lTohlinell2002h . These imply 
a smaller X-ray emissivity per unit stellar mass by 
CVs and ABs, rather than a higher emissivity. 

3.2. Magnetic Field Configuration 

The most puzzling aspects of the GCDX is 
its high temperature. Since the kT w 7keV 
plasma is too hot to be gravitationally bound, it 
requires a huge energy source without the con- 
finement of the plasma. One idea to address 
this energetics problem is the c onfinement of the 
plasma by magnetic fields (MFs; iMakishiina 1994t 
Tanuma et al. 1999() . If a large-scale toroidal 
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MF is developed and sustained, and the MF is 
strong enough for nearly energy equipartition 
with the plasma, the GCDX could be almost con- 
fined within the NB. However, the large scale MF 
configuration was thought to be predominantly 
vertical, suggesting that the magnetic confine- 
ment does not work well, although observations 
of the MF in the NB have been limited to the 



regio n in d e nse molecular clo uds (jNovak et al 



2OOOI l2003t IChuss et all l2003l) and verv thin 
non-thermal radio filaments (ITsuboi et alj Il986t 



Yusef-Zadeh et al.lll997t iLang et al lll999[ ). Recent 
observations have revealed that NIR and wide- 
field polarimetry offers a promising tool to trace 
a large-scale MF, and a troidal configuration near 
the Galactic plane has been claimed in the GC 
(jNishivama et al.ll2010l) . 
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Fig. 5. — Polarimetry results covering 3?0 x 2?0 
in the Galactic coordinate, toget her with an inten- 



sity m ap of 6.7 keV line emission ([Nobukawa et al 



20121 ). The cyan vectors show the inferred mag- 
netic field direction, and the lengths are propor- 
tional to polarization percentage. The vectors are 
averaged in a circle of 2'A radius with a 3.'0 grid, 
and plotted with thick bars (detected with more 
than 3(t) and thin bars (detected with 2 — 3ct). 

The obtained polarization map (Fig. [S]) sug- 
gests a large-scale toroidal MF configuration in the 
NB. The histogram (Fig. ^ of the MF directions 
at I 6 |< 0?4 has a clear peak at 90° which is the 
direction parallel to the Galactic plane. On the 
other hand, at high Galactic latitude {\b\ > 0?4), 
the fields are nearly perpendicular to the plane, 
i.e., poloidal configuration. This suggests a tran- 
sition of the large-scale configuration, and such 



a transition can be naturally explained by the 
time evolution of MFs. An initially predominantly 
poloidal, larger-scale MF is sheared out in the az- 
imuthal direction b y the differential ro tation of an 
accreting gas disk (|Uchida et al.lll985() . The tran- 
sition region, b ~ 0?3 — 0?4, is in good agreement 
with the scale height o f the 6.7 keV emission, 0?27 
(jUchivama et al.1120131 ). 
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Fig. 6. — Histograms of the magnetic field direc- 
tion at I 6 |< 0?4 (red), | 6 |> 0?4 (blue), and both 
(black). The red histogram has a clear peak at 
the direction almost parallel to the Galactic plane 
(90°), while the blue one has a peak at ~ 150°, 
almost perpendicular to the plane. 

The determination of the MF strength is still 
quite difficult in this region, but it seems to 
conve rge to the val u e of 50 < BjuG] < 200 
(e.g., iFerrierd bOlll ICrocker et al] boill) . At 
this field strength, the magnetic energy density 
(~ 0.1 — IkeVcm"^) reaches nearly equipartition 
with those of diffuse hot plasma (~ 0.5 keV cm~'^) 
and gas turbulence (see Fig. 4 in ICrocker et al 



20101) . This suggests that MFs provide significant 
pressure support against the diffusion of the hot 
plasma. 

If the plasma were not supported, it would 
be rushing out of the Galactic plane vertically 
as a galactic wind. The escape velocity, typ- 
ically several hundred kms~^, is smaller than 
the sound speed of the 7keV hot plasma of 
1, 400kms~"'^. Assuming the gas flows out from 
the X-ray emitting region at the sound speed, the 
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escap e timescale is 4 x 10'' yr (jBelmont et al 
20051 ) ■ This requires a huge energy input to sus- 
tain the hot plasma; e.g., an unre ason ably high su 



perno va rate of ~ 5 x 10 yr ^ (|Uchivama et al 



20131 but note that kT « 7keV plasma is hot- 
ter than is observed in supernova remnants). If 
the plasma is magnetically confined, and there is 
no other cooling mechanism, the hot plasma only 
cools by radiation with a timescale of 10^ — 10* yr 
(jMuno et al.ll2004l) . several orders of magnitude 
longer than the escape timescale. This would re- 
duce the required energy input by several orders of 
magnitude and thus relax the energetics problem. 

There is no widely accepted mechanism to heat 
the plasma to kT ss 7keV. If past activities o f 
the supermassive black hole (jKovama et al.lll996l ) 
heat ambient interstellar gas, created plasma 
could be confined magnetically. If a troidal MF 
is sustained by a differential rotation of a gaseous 
disk, magnet ic reconnection is a possible heating 
mechanism ( Tanuma et al. Il999l) . Star forma- 



tion activities and resulta nt supernovae are also 
implied to be the origin (jCrockeii 120121 ) , while 
supernova-driven outflows advect materials out of 
the GC. However, outflows from the GC seems 



to be magnetized ijCarretti et al.ll2013l ). hence the 
existence of a large-scale troidal field and outfiows 
might not be contradictory. 

4. Summary 

We have used imaging and polarimetric data 
sets of the GC region to investigate the origin of 
the GCDX. We have constructed a stellar number 
density map, and compared its longitudinal and 
latitudinal profiles with those of the 6.7keV emis- 
sion. We have estimated that the contributions 
of the old stellar population to the GCDX at the 
NSD and NSC are - 50% and - 20%, respec- 
tively. Our findings support the notion that the 
GCDX is not only caused by a population of un- 
resolved point sources but must also stem from a 
hot interstellar plasma component. Polarimetric 
observations reveal a large scale toroidal magnetic 
field configuration which allows a magnetic con- 
finement of the hot plasma. 
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